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The diffusion of hard-core particles subject to a global bias is described by a nonlinear, anisotropic 
generalization of the diffusion equation with conserved, local noise. Using renormalization group 
techniques, we analyze the effect of an additional noise term, with spatially long-ranged correla- 
tions, on the long-time, long-wavelength behavior of this model. Above an upper critical dimension 
dLR, the long- ranged noise is always relevant. In contrast, for d < dtjj, we find a "weak noise" 
regime dominated by short-range noise. As the range of the noise correlations increases, an intricate 
sequence of stability exchanges between different fixed points of the renormalization group occurs. 
Both smooth and discontinuous crossovers between the associated universality classes are observed, 
reflected in the scaling exponents. We discuss the necessary techniques in some detail since they are 
applicable to a much wider range of problems. 
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I. INTRODUCTION 



Systems coupled to multiple energy reservoirs, sustaining net transport currents, are prevalent in nature but fall 
outside the fundamentally well-understood paradigms of equilibrium statistical mechanics. The characterization and 
classification of such systems, starting from simple microscopic models, remains a key problem of current research. 
Since steady states in such systems constitute the closest relatives of equilibrium (Gibbs) states, their study has 
attracted much recent attention [Q, revealing a wealth of unexpected phenomena, even in simple model systems. One 
of the prime, and most elementary, examples of this kind is a system of hard-core particles undergoing biased diffusion. 
With appropriate boundary conditions, this system, also known as the simple asymmetric exclusion process (ASEP) 
[2 1, can exhibit a number of surprising features, including anomalous diffusion long-range spatial correlations 

[5|, and shocks ||. One of its key characteristics is the breaking of the usual fluctuation-dissipation theorem [Q. 

In the following, we will discuss the effect of long-range correlated noise on this model, formulated via a Langevin 
equation in continuous space and time. Generated, for example, by correlated fluctuations in the strength of the 
bias, the additional operator acts in the one-dimensional subspace (labelled "longitudinal") selected by the bias. Its 
momentum dependence takes the form gii 2 *- 1- "', with < a < 1. The two extreme values, a = and a = 1, 
are of particular interest. The choice a = corresponds to the standard ASEP. In contrast, a = 1 leads to a 
Langevin equation with non-conserved noise, describing biased diffusion of particles which can occasionally be created 
or annihilated. First introduced by Hwa and Karclar jq] as a continuum description for sliding avalanches in sandpilc 
models, it was analyzed in more detail by Becker and Janssen |)|. 

The key advantage of the correlated noise term is that it allows us to capture the crossover between these two 
models, characterized, respectively, by conserved and non-conserved noise. Thus, using field theoretic techniques, we 
discuss the full renormalization group (RG) flow in (a, d) space, for < a < 1. While the limit a — > 1 presents no 
difficulties, the opposite limit a — > is far more subtle: we will see that, just below the upper critical dimension 
d c (a), two nontrivial renormalizations are required to render the a = theory finite, in contrast to a single one if a 
is positive and of O(l), leading to an apparent discontinuity in the critical exponents. A similar situation arises in 
standard </> 4 -theory with long-ranged interactions JlO| . There, however, the discontinuity is entirely spurious and can 
be removed: letting e = d c (0) — d denote the distance from the upper critical dimension of the short-range theory, 
the key is to recognize |T^] that there is a region of small a = O (e) , where a careful analysis of the RG flow reveals 
a smooth crossover between the a = and the a = O(l) theories. Thus, all universal scaling properties are shown 
to depend continuously on a and d. Here, in contrast, the limit a — > is even more intricate: while the short-range 
theory is controlled by two fixed points, only a single one remains in the long-range model. Thus, only one of the two 
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short-range fixed points is smoothly connected to the long-range fixed point, leading to a continuous crossover. The 
other short-range fixed point eventually loses its stability, accompanied by a discontinuity in critical exponents. 

A second feature of our model concerns its close relationship to the KPZ |fl| or Burgers' |12) equation with correlated 
noise [ p^|Jl4| . It is well known that the continuum description for the ASEP and the noisy Burgers equation coincide 
in one spatial dimension. In contrast to the Burgers equation, the £-expansion for biased diffusion around d c = 2 
presents no difficulties, and numerous results can be obtained exactly, to all orders in perturbation theory. These may 
therefore be analytically continued to d = 1 and are expected to hold for both, the ASEP and the KPZ equation. In 
fact, to the extent that exact solutions flq| are available, this expectation is indeed confirmed. Thus, our analysis 
bears some relation to the behavior of the one-dimensional Burgers equation with correlated noise. 

This paper is organized as follows. We first introduce the Langevin equation for the "short-range" (SR) version 
of our model, corresponding to the usual c£-dimensional ASEP with translational invariance. We then add a second 
noise term with long-range spatial correlations and propose two possible microscopic mechanisms for such a term. 
Turning to the RG analysis, we first summarize the SR case, a = 0. The "long-range" (LR) theory associated with 
a = O (1) is presented next, and the limit a — > 1 is discussed. Finally, we introduce the "hybrid" model, characterized 
by a — O (e), which interpolates between the SR (a = 0) and LR (a = 0(1)) theories. Its RG flow is computed in 
a double expansion, where both a and e are small parameters. The different fixed points are interpreted and their 
stability is evaluated, illustrating how the subtle crossover from a = to a = 0(1) occurs. We conclude with some 
comments and open questions. 



II. THE MODEL 



We consider a <i-dimensional system of hard-core particles, which are allowed to diffuse on a regular lattice with fully 
periodic boundary conditions. The bias, reminiscent of an electric field E = Ee^ acting on charges, favors particle 
moves along a specific ("longitudinal") direction along unit vector e||. The long-time, long-wavelength properties of 
this model are most conveniently captured by a coarse-grained description in continuous space and time. Here, only 
a single slow variable is needed, namely, the local particle density c(r, t). Since the associated Langevin equation 
has been discussed previously [HJJ], we describe it only briefly here. Due to particle number conservation, it takes 
the form of a continuity equation, dtc + Vj = 0. In addition to a diffusive term, the current j contains an Ohmic 
contribution, j_E = k (c) E, induced by the bias, and a Gaussian noise term, modelling the fast degrees of freedom. By 
virtue of the excluded volume constraint, the conductivity n (c) must vanish with both particle and hole density, i.e., 
k (c) oc c (1 — c). Aiming for a perturbation expansion around a field with zero average, we introduce the deviations 
s (r, t) = c (r, t) — c from the average density c. A term linear in s, arising from the expansion of n (c), can be absorbed 
via a Galilei transformation. Finally, we allow for anisotropics in the diffusion coefficient and the correlations of the 
Langevin noise, since the bias singles out a preferred direction. Neglecting contributions that are irrelevant in the 
long- wavelength, long-time limit, the Langevin equation can be summarized in the form Q: 

d t s = A { (Vi + P d 2 ) s + |as 2 } + r, (r, t) . (1) 

Here, Vj_ (d) denote the gradients in the transverse (longitudinal) subspace, and the kinetic coefficient A defines a 
time scale. The coupling g oc E captures the effects of the drive. The noise 77 (r, t) is Gaussian so that two moments 
suffice to characterize its full distribution (after a suitable rescaling of the variables) : 

(V(r,t)) =0, 

(r, (r, t) rj (r', f)> = -2A (Vi + ad 2 ) 5 (r - r') 6 (t - f ) . (2) 

The differential operator (V^_ + crd 2 ) in the second moment ensures that the conservation law is strictly obeyed. 
Clearly, the correlations described by this form are purely local, so that we will refer to Eqns. ([l]J|) as the "short- 
range" (SR) theory. Its universal properties have been discussed in ||,|5| . 

The simplest way of introducing a correlated noise into this equation is to add a long-range correlated noise to the 
local one. Thus, Eqn. @ is amended to 

(V(r,t)) =0, 

(tj (r, t) 77 (r, t')) = 2A [- (Vi + ad 2 ) + b (-0 2 ) '""j 5 (r - r') 5 (t - t') . (3) 

The new operator b (— d 2 ) 1 a reflects a power-law decay of the noise correlations in real space, giving rise to a q 21 " 1 °^ 
momentum dependence in Fourier space. Clearly, setting a to zero reproduces the SR theory, albeit with a replaced 
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by a + b. In contrast, choosing a = 1 generates a non-conserved noise, playing the role of a random source term in 
Eqn. (|l|). Note that the vanishing of the first moment ensures that the particle density, while not strictly invariant 
under the dynamics, remains conserved on average. The long-time, long- wavelength behavior of this theory has been 
analyzed in ^M. Other choices of long-ranged noise correlations are of course possible |1(J. The version adopted 
above is "minimal" in the sense that, first, it leads to a nontrivial, renormalizable theory, and second, it allows us to 
discuss the full crossover between the conserved (a = 0) and non-conserved (a = 1) cases using just one additional 
operator. Microscopically, such long-ranged contributions to the Langevin noise can be generated in two ways: one 
option is to impose a bias which consists of a spatially uniform and a locally random component, with the latter 
controlled by a distribution with long-range correlations. Alternatively, one can retain the strictly uniform bias but 
randomly add or remove particles, according to a long-ranged distribution. More details will be presented in [|l6). 

Due to the presence of the nonlinearity, the computation of averages over the noise distribution requires a pertur- 
bative approach, coupled with renormalization group methods. For these purposes, it is most convenient to recast 
the Langevin equations as dynamic functionals p7| . Introducing a Martin-Siggia-Rose |18j response field s(r,t), we 
can recast Eqns. (0J|) in the form 



J [8,8] = I dt I d d r{§d t s-X§ (Vl + pd 2 ) s + |<9s 2 



+XS 



ivl + ad^-bi-d 2 ) 1 - 



(4) 



so that correlation and response functions can be computed as functional averages with weight exp (-J). We stress 
that the SR model is recovered from this general expression, simply by setting 6 = 0. 

As a first step towards the RG analysis of this model, we discuss its scaling symmetries. First, under a global 
rescaling of all coordinates, t^— > /z r*, we find a scale invariant theory (as long as we neglect any cutoff-dependence) 



provided Xt 



2 Xt, s -► n<.*-**)/* a an d § 



M (d+2a)/2~_ 



Moreover, we obtain g — > fx 



(l+a)-d/2 



g which allows us 



to identify the upper critical dimension d c (a) = 2 (1 + a). To avoid confusion, we will reserve the symbol e for the 
SR case, i.e., e = 2 — d. For the general case, we define e = d c (a) — d — 2 (1 + a) — d. Second, we may rescale 
the parallel coordinate alone, due to the anisotropy of the model: under x\\ — > f3x\\, scale invariance is retained if 



s -> /T 1/2 s, 3 f3~ 1/2 s, p 
of effective couplings, namely, 



u = G d g 2 p 



f3 2 a, b — > /3 2 ( 1_Q ^6 and g — > P 3 ^ 2 g. This permits us to identify a set 
~ 5 / 2 <7, w — pja and v = bp a /a. These will emerge quite naturally in 
perturbation theory, combined with appropriate powers of fi which absorb the remaining momentum dependence. 
For convenience, a geometric factor Gd = Sa/ (27r) d , with Sa the surface of the d-dimensional unit sphere, has been 
absorbed into u. For later reference, we note that only the positive octant u > 0, w > and v > corresponds 
to physical theories even though the RG flow can be discussed in a larger parameter space. Finally, the dynamic 
functional (||) exhibits a nontrivial continuous symmetry, with parameter a, which is characterized by the "Galilean" 
transformation s (r, t) — ► s (r + Xage^t, t) + a, s (r, t) — ► s (r + Xagent, t) . 

So far, we were able to consider the case of general < a < 1. The next section will show, however, that this 
procedure cannot be continued for the actual perturbative calculation of correlation and response functions. 



III. RENORMALIZATION GROUP ANALYSIS 



We now turn to perturbation theory, with the goal of identifying those correlation and response functions that 
require renormalization, due to the presence of ultraviolet (UV) divergences in Feynman diagrams. As usual, we focus 
on the one-particle irreducible (1PI) vertex functions with n (n) external s- (s-) legs, Tn n . These will be computed 
in dimensional regularization using minimal subtraction. Since the models associated with a = and a = 1 have 
been discussed previously (in and ||,|9|], respectively), we only briefly summarize their RG structure, focusing 
predominantly on the theory with < a < 1. 



A. The short-range theory: a = 0. 

We first review the model with a — PP]. Here, the coupling b and the effective v are redundant and will be set 
to zero. The upper critical dimension is d c (0) = 2. Straightforward dimensional analysis shows that there are three 
naively divergent vertex functions: Tn, and T^. A Ward identity based on the "Galilean" symmetry shows, 
however, that T12 requires no additional renormalization once Tn and have been rendered finite. Letting the 
superscript ° denote bare couplings, the renormalized couplings are defined in the usual way: p — Z^ R p, a — Zp R <r, 
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g = Z g g(i £ / 2 , where e = 2 — d. The parameter p, sets a typical momentum scale which will control the RG flow. All 
SR functions are marked with an explicit superscript to distinguish them from their finite a counterparts. The UV 
divergences in Tn and I^o give rise to nontrivial Z-factors for p and a, which can be computed perturbatively in u. 
A one-loop calculation yields Z s p R = u(3w + l)/ (16e) + O (u 2 ) and Zf R = u (3w 2 + 2w + 3) / (32e) + O (u 2 ). We 
emphasize that w is not treated perturbatively here; therefore, O (1) fixed points for w should not come as a surprise 
later. Finally, Z g = 1 to all orders in u, due to the "Galilean" invariance. To complete the discussion, we introduce 
the Wilson functions 



^ R =pd„hia\ 



3w + 1 , 9 > 
-«^6~+0(« 2 ). 

3w 2 + 2w + 3 n , 2 v 
-U + 0(u) 



(5) 



The RG flow is expressed through the ^-dependence of the renormalized couplings u = Gag and w = p/o~, 

especially in the scaling limit p, — > 0: 



bare 



3(7 



(0) 



In this form, the right hand sides of Eqn. (g) are exact to all orders. 

These flow equations possess three fixed points, marked by the vanishing of both f3^R and (3^ R : (i) u* = 8e/3 + 
O (e 2 ), w* = I; (ii) u* = 16e + O (e 2 ), w* = 0; (hi) m* = 16e/3 + O (e 2 ), w* = 1/3 + O (e); and finally a fixed line 
(iv) u* = 0, with arbitrary w. The stability of these fixed points (and line) can be expressed through the 2x2 matrix 
M-Sr = (di0j RS ) , i,j = u,v, of derivatives of the /3-functions, evaluated at the fixed points. We find that the fixed 
line (iv) is stable only for e < 0, i.e., d > 2, and will therefore be labelled a "Gaussian" line. Focusing on e > 0, we 
find that (i) and (ii) are both locally stable. Fixed point (iii) is hyperbolic and sits on the separatrix w = 1/3 + O (e) 
which separates the domains of attraction of (i) and (ii). Flow lines starting in the region w > 1/3 + O (e) are 
attracted towards fixed point (i), and vice versa. We refer to fixed point (i) as the FDT- restoring short-range fixed 
point, since theories with w = 1 can be shown to possess a higher degree of symmetry, associated with satisfying 
the FDT Q Here Eqn ••'!'•<<•- ><- i.. ..,.•....»!.)<. ii>„ Hv,»i -,,]-,.,. ,,r . sii ,„„i , sh 

(SR* = (SR* _ _2 £ /3. Conversely, fixed point (ii) violates the FDT, and (jj 
are known only perturbatively. 

Finally, we comment briefly on the scaling properties associated with these fixed points pLS], which arise from 
the RG equations for the vertex functions. For example, the dynamic density correlation function (structure factor) 
S (q, uj) = (s (— q, — lo) s (q, w)) scales as 



allows us to compute the fixed point values of (p and (% R exactly, to all orders in e: 



»SR* 



e + 0(e 2 );(S 



SR* 



-3e/2 + Q(e 2 



-1-A. 



«, q± r\u;r 2 ) 



(7) 



This reflects the emergence of anomalous diffusion, characterized, near fixed point (i), by a strong anisotropy exponent 
§ A S r = -(% R */2 = (2 - d) /3, and, near fixed point (ii) by an exponent A S r = -^ SR * 
we turn to the case of finite a. 



(SR*/ 2 = e/2 + 0(e 2 ). Next, 



B. The long-range models with a = O(l). 

As indicated in the introduction, we will have to distinguish between values of a which are O (2 — d) versus those 
which are O(f). Here, we analyze the second case, which corresponds to true long-range (LR) theories. The first 
case, being a "hybrid" between short-range and long-range models, will be deferred to the next subsection. 

Recalling our discussion in Section 2, the critical dimension is now d c (a) = 2 (1 + a), and we define e = 2 (1 + a) — d, 

to be distinguished from e = 2 — d. The operator s (Vj_ + ad 2 ) s is clearly irrelevant compared to bs (—d 2 ) " s and 
may be dropped. The naive dimensions of the fields s and s arc a-dependent, and as a result only two vertex 
functions, namely Tn and Ti2, are naively divergent. In contrast to the short-range case, is naively convergent 
here. Moreover, all divergent contributions to any vertex function are polynomial in the momenta [[l0|| , so that b is 
not renormalized. The "Galilean" invariance still holds, so that Ti2, and hence g, requires no new Z-factor. Thus, 
defining renormalized couplings according to p = Z^ R : b = Zt,b, g = Z g gp £ / 2 , we find Zb = Z g = 1 to all orders, and 
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Zp R = 1 — uA (a) je + 0(u 2 ), to first order in u = uv. The latter is the appropriate effective coupling here, since a 
no longer appears in the theory. The coefficient A (a) = T (1 + a) T (3/2 — a) (1 + 2a) j (8->/W) controls the one-loop 
divergence in Tn. The Wilson functions are easily obtained: 

Q L p R ee [idfj_ lnp\ bare = ~uA (a) + O (u 2 ) (8) 

and 

'5 



fc«EE^u\ bare 



LR 



(9) 



Clearly, the LR theory exhibits one (Gaussian) fixed point u* — which is easily shown to be stable above the upper 
critical dimension, and a nontrivial fixed point u* — 2ej [(5 — 2a) A (a)] +0 (e 2 ), stable below d c . At this fixed point, 
the value of is again known to all orders in e, namely, ( p R * = — 2e/ (5 — 2a). 

The scaling properties of the long-range models are again easily derived from an RG equation for the vertex 
functions. For comparison with the SR case, we quote the result for the dynamic structure factor: 

S (q, u) = r 2 ^s (giii- 1 -***, qi_l~\ u;r 2 ) . (10) 

Once again, we observe anomalous diffusion; here, however, it is controlled by the LR exponent A^r ee —( p R */2 = 
(7-d)/(5-2a)-l. 

It is obvious at this point that the limit a — ► 0, taken in the Wilson function (||) or the correlation function (|lC|), 
will not restore the SR theory. First, lim Q ^o Cp R = —u/16 + O (u 2 ) which reproduces only the w = component of 
Cp R - Second, this naive limit cannot generate a non-vanishing £^ R which plays a key role in the SR theory. Thus, not 
surprisingly, we do not observe a smooth crossover from the SR to the LR theory, by simply letting a tend to zero 
in this naive fashion. An elegant way of dealing with this difficulty was suggested in |l(J: By treating a as a small 
parameter of O (e) and expanding in both e and a, one can "magnify" the small a-section of (d, a) space and resolve 
the crossover between the SR and LR fixed points. Concerning the scaling form of the structure factor, the prefactor 
l~2(i+a) or igj na tes entirely in the bare dimensions of the fields. Since these are well defined under a naive a — > 
limit (cf. Section 2) and remain unrenormalized, we have lim^^o ^~ 2 ' 1+a ' — l~ 2 without further complications. The 
exponent Alr, however, must be tracked much more carefully. 

Before turning to these subtleties, we briefly consider the other limit, a — ► 1. First, we summarize the a — 1 theory 
pL with e = 4 — d, in its own right. Here, b measures the strength of the (non-conserved) noise. The effective coupling 
is u — Gdg 2 p~ 3 ^ 2 b. According to there is only one nontrivial Z-factor, namely Z p = 1 — 3m/ (8e) + O (u 2 ). Again, 
Zb = Z u = 1 to all orders. Thus, the theories with a = 1 and < a < 1 exhibit the same set of UV divergences. 
From the preceding discussion, it is therefore clear that we may anticipate a smooth crossover here. This is explicitly 
confirmed by the results of || for a — 1 : two fixed points u* = and u* = 16E/9 + (e 2 ) are found, stable above and 
below d c = 4, respectively. At the nontrivial fix ed p oint, A = — Cp/2 takes the exact value e/3. It is easy to check 
that this agrees with the a — > 1 limit of Eqns. (p||l0|), so that this limit is continuous. 



C. The hybrid theory: a = O (e). 

Finally, we turn to the analysis of the key region in (d,a) space, namely, a = 0(e). The naive a — > limit 
presupposes £ a, and hence fails to resolve the crossover between the SR and LR theories which occurs for 
a = 0(e). With both a and e small, we follow |l(]] and analyze the general dynamic functional, Eqn. (ji|), near the 
upper critical dimension of the SR theory, i.e., d = 2. Both couplings a and b will be retained, and their RG flow 
will be studied. We refer to this model as the "hybrid" theory, since it contains the vital elements of both SR and 
LR cases. Thus, the hybrid theory will exhibit a well-defined a — > limit. Once again, we will be able to obtain a 
number of results to all orders in e — 2 — d. 

Our first task will be to identify a set of suitable couplings which allow us to interpolate between the SR and the 
LR theories. Guided by JlOj ] , we consider the general structure of the Wilson functions in the SR and the hybrid 
model. The appropriate effective couplings will be those that map those functions onto one another. 

To obtain a general form for the Wilson functions, we begin with the Z-factors and construct them order by order 
in perturbation theory. Focusing on the SR theory first, a typical graph of Tn, at L-loop, contains 2L vertices, L 
(bare) correlators and 2L — 1 (bare) propagators. Each vertex carries a factor of gq\\ , and each correlator contributes 

a factor (^qj^ + aq 2 ^, where (u>; q) = (cj;q'||,qx) denote the frequency and momentum of the corresponding line. 
Symmetry requires that 2, out of the L parallel momentum factors contributed by the vertices, will not be integrated 
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over. Moreover, all denominators are generically of the form (j,ujX 1 + q\ + PSfV Changing the parallel integration 
variable from q\\ to p 1 ^ 2 q\\ , it is straightforward to show that such a graph will give rise to a contribution of the form 

u L J2m=o wm ALm/ {Le) + O (e~ 2 ) in Zp R . Here, the A^ m are the numerical coefficients of the O poles, in 

minimal subtraction. Higher order poles appear also, of course, but cancel in the Wilson functions and need not be 
considered for this reason. Generic i-loop graphs of I^o have a similar structure, with L + 1 correlators and 2L — 2 
propagators, and the coefficients of the lowest order e-poles will be denoted by -B_Lm here. Collecting, we obtain the 
Z-factors in the general form: 

oo L . 

^ = i+E^E- m %+°(^ 2 )> 

L=l m=0 

Z™ = l + J2u L ^w™^+0 (e- 2 ) . (11) 

L=l m=0 

The general form of the Wilson functions is easily found from the above: 

oo L 

L—\ m=0 
oo L + l 

C S a R = E u L £ w m B Lm . (12) 

L = l ro=0 

Next, we consider the hybrid case. Two key differences emerge: First, each correlator now contributes a factor of 
a given diagram. This product gives rise to a sum of individual terms each of which contains 

L — n factors of type [q\ + crq^j and n factors of type bq 2 ^~ a \ with n = 0, 1, L. Second, since a is now of O (e), 

the poles in these individual contributions are proportional to 1/ {Le + 2na) [H . Following the same reasoning as 
above, and recalling the effective coupling v — bp a /a, one finds quite readily that 

oo L L — n / T \ . 



p l^l^\ n ) Le + 2na 

L = l n=0m=0 v 7 
oo L+l L+l — n ,j s „ 



L = l n=0 m=0 



Here, A^ m „ and Bhmn are the numerical coefficients of each pole. In minimal subtraction, they are independent of 
both e and a, since the latter is also 0(e) here. Clearly, the associated Wilson functions take the form 

oo L L—n / t- \ 

L=l n=0ro=0 



L + l L+l-n 



E E (^V^- (i4) 



E<<'- v " ^ 

L=l n=0 m=0 



Next, we establish a relationship between the coefficients (Ai, m , Bi m ) and (Ai, mn ,Bi, mn ), by considering the a — > 
limit of the hybrid theory. In this limit, the dynamic functional, Eqn. (0), simply reduces to that of the SR theory, 
with a replaced by a + b. Therefore, we can relate limc^o Z p and lim^^o Z a to Z p R and Z^ R . Due to the replacement 
of <t by (7 + 6, the SR Z-factors depend on a modified set of couplings. Defining 

U = u (1 + v) , w 



1 + v' 
we obtain: 

lim Z p (u, w, v) = Zf R (u, w) , 

a — >0 



lim Z a (u, to, v) = (1 + v) Z s a R (u, to) - v. (15) 
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Inserting the explicit forms for the Z-factors into Eqn. ( |15[ ) and recalling that (Ar, OT „, Br, mn ) are independent of a, 
we read off the relation between the two sets of coefficients: 



L 

n 

L + i 

n 



L — m 
n 

L + 1 — m 

n 



Bhr, 



(16) 



Combining Eqns. ( |l6| ) and ( [l4| ) provides us with an identity between SR and hybrid Wilson functions: 

C p (u,w,v) = Cp R (u,w), 

Cr(u,W,v) = (1 + V){S R (U,W). 



(17) 



Here, (p R and ^ R denote the Wilson functions of the SR model, Eqn. (||), evaluated at (u,w). In the following, we 
use the abbreviated notation £ p = Cp R (u, w) and = Cr R (u, w). 



Eqns. (|17]) form the basis for the remainder of the paper: they hold the key for the discussion of the hybrid theory 
and for the desired matching between LR and SR models. It is particularly gratifying that they are valid to all orders 
in perturbation theory. An explicit calculation of the hybrid Wilson functions, instead of the general considerations 
presented above, would establish Eqns. Jl7|) only up to a given order. 

To discuss the RG flow for the hybrid theory, we compute the 0- functions for the couplings (u,w,v): 



f3 a (u,w,v) = ^u\ bare 
f3yj (u,w,v) = ndpw\ bare 

f3 v (u,w,v) = ndf,v\ bare 



2-Cp 



Cp - Cx + ol- 



(2-Cp) 



1 + v 

(l+v)( a + a(2-t p ) 



UK 



(18) 



and seek their fixed points (u* , w* ,v*). These fall into two groups. 

The first group is characterized by v* = 0: In this case, Eqns. ( filf ) reduce to the fixed point equations @ for 
the SR theory, so that we recover the equivalents of the familiar SR fixed points, namely (i) the FDT satisfying 
(8e/3 + O (e 2 ) ,1,0), (ii) the FDT violating (l6e + O (e 2 ) ,0,0), (hi) the fixed point (16e/3 + O (e 2 ) ,1/3 + (e) , 0) 
on the separatrix between (i) and (ii), and finally the Gaussian fixed line (iv) (0,w, 0). In contrast to the SR case, 
however, their stability needs to be investigated in the i/iree-dimensional space spanned by (u, w, v), controlled by the 
3x3 matrix M = (di(3j)* , i, j = u,w, v. For v* = 0, we find dufi v = d^Pv = at all of these fixed points, so that only 



the upper left 2x2 corner of M (which is just 



and the bottom right element, d v /3 v = — a (2 — £ p ), are 



relevant. The eigenvalues of M sfl were already computed in Section (3.1), so that only fixed points (i) and (ii) remain 
as candidates for global stability if e > 0. For fixed point (i), Q — Q — ~2e/3 to all orders (cf. Section 3.1) so that 
d v (i v is positive for a < a.\ = ej (3 + e). Thus, a\ demarcates the stability boundary of the FDT satisfying SR fixed 
point (i). Considering fixed point (ii), we have Q — — e + O (e 2 ) and Q* = — 3e/2 + O (e 2 ), giving positive values for 
d v (3 v provided a < a 2 = 3e/4 + (e 2 ). So, the FDT violating SR fixed point (ii) remains stable for a larger region of 
a than its FDT-satisfying partner. We remark that, in contrast to ai, a 2 is known only perturbatively. Finally, the 
fixed line (0, w, 0) is stable provided both e and a are negative. 

Returning to e > 0, it is clear that another fixed point must become stable beyond a 2 - By necessity, this can only 
be a member of the second group, having v* ^= 0. For all of these, the [. . .] bracket in the last line of Eqns. ( flg| ) 
vanishes, so that the first two equations simplify to 



Pu (u,W,V*) 



Pu, (U,W,V*) = 



(e + 2a)+[--a)C P 



2a + (1 - a) Cp 



u. 



(19) 



One should note that, with e — e + 2a, (3u is precisely the /3-function of the LR theory, Eqn. (j^). Thus, we identify the 
second group as the hybrid partners of the LR fixed points. The equation for p a has a trivial solution u* = and a 
nontrivial one, with u* ^ 0. Seeking the corresponding values for w* , we obtain a "Gaussian" fixed point (v) (0, 0, oo) 
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and a nontrivial one (vi) with u* = 32 (e + 2a) /5 +0 (e 2 ) , w* = and v* = (4a - 3e) / [3 (e + 2a)] + 2a/3 + (e 2 ) . 
Note that we have expanded u* in both s and a, keeping only terms to first order in cither, since both are assumed to 
be of the same order, v* is 0(1) which should not be disturbing since v is not a perturbative coupling. Since u* ^ 
here, we find (* = -2 (e + 2a) / (5 - 2a) to all orders. 

Turning to the stability of these fixed points, it is straightforward to determine that the Gaussian line is stable for 
d > 2 (1 + a), i.e., above the upper critical dimension of the LR theory. For d < 2 (1 + a), we find a more complex 
situation: for a < a\, fixed point (vi) has two unstable directions. One of these becomes stable above a\. For a > at2, 
this fixed point is globally stable. 

It is natural, of course, to seek fixed points where neither u* nor w* vanish. This attempt gives us a pair of 
equalities, valid to all orders: 

0=(e + 2a)+Q-a)c;, 

= 2a+(l-a)C;, (20) 

which do not result in a unique equation for such a fixed point. Instead, they select a specific 
value of a, namely a — a\ = e/(3 + e), where a fixed line (vii) exists, parameterized by w: 
(32e/ [3 (3w + 1)] , w, (3w — 1) (1 — w) / (3m) 2 + 2w + 3)) . This line mediates the stability loss of fixed point (i). 

The crossover scenario between the SR and the true LR theory can now be summarized. Let us fix d just below 2 
and increase a starting from zero. For sufficiently small < a < a\, the RG flow is dominated by the SR behavior. 
The SR fixed points (i) and (ii) are stable, within their respective domains of attraction. The separatrix forms a 
surface which cuts the v = plane at w = 1/3 + O (s) and then bends over to larger values of w as v increases. The 
unstable LR fixed point (vi) is found on the separatrix, in the unphysical region v < 0. When a reaches a±, the 
correlated noise begins to make its presence felt. Specifically, for a\ < a < a 2 , the FDT violating fixed point (ii) 
is the only globally stable fixed point. The FDT restoring fixed point (i), while still stable within the v = plane, 
has become unstable to small perturbations out of that plane. Thus, a flow line starting near (i), with a small v > 
component, will first flow out into the half-space v > and then bend back towards v — 0, flowing into the FDT 
violating fixed point (ii). The sign of v remains invariant under the flow. The LR fixed point (vi) is still unphysical, 
but has moved closer to the v — plane. Finally, at a — ct2, (vi) merges with (ii) and moves out into the positive 
v region as a increases beyond oli- The LR fixed point (vi) is now the only stable one, and the global RG flow is 
dominated by the LR theory. A different view of the SR-LR crossover is presented in Fig. 1 which shows the location 
of the stability boundaries as functions of a and d. 




1/4 1/2 3/4 1 

a 



FIG. 1. Stability boundaries of different fixed points as functions of d and a. The heavy solid lines denote boundaries whose 
location is known to all orders in perturbation theory. 
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It is interesting to note the two different scenarios which control the stability loss of fixed points here. Fixed points 
(ii) and (vi) exchange their stability by merging, similar to the stability exchange between the Gaussian and the 
Wilson-Fisher ^| fixed points in ^-theory when d drops below 4. In contrast, fixed point (i) loses its stability, at 
a = oti, by generating a fixed line. This fixed line has two key properties: first, it connects fixed points (i), (iii) and 
(vi); second, it lies on the surface forming the separatrix between fixed points (i) and (ii). This surface bends back 
towards the v — plane, as a increases, to the extent that it reaches fixed point (i) at a = a%. Thus, fixed point (i) 
loses one of its three stable directions when being absorbed by the separatrix; moreover, this direction is now spanned 
by the fixed line: this allows the flow to change sign when a crosses ax- 

The crossover which we observe in the flow diagrams is also reflected by the exponents. As an example, we consider 
the scaling form of the LR structure factor, Eqn. (|l0|). We recall that its scaling behavior is determined by the 
exponent A^r = s/ (5 — 2a). The hybrid fixed point (vi), on the other hand, generates an exponent A = — £*/2 = 
(e + 2a) I (5 — 2a). These two exponents clearly match, since e — e + 2a. As a decreases, we reach the stability 
boundary a 2 = 3e/4 + 0(e 2 ), where (vi) and (ii) exchange stability. Here, the exponent associated with (ii) is 
Asr — —(p R */2 = e/2 + O (e 2 ), to be compared with the corresponding value near (vi), A = (e + 2a 2 )/(5 — 2a 2 ). 
In this way, continuity is ensured. However, a discontinuous change of exponents may occur as a decreases below 
a\. above ax, all positive (u, w, v) fall into the domain of attraction of fixed point (ii), with Asr = e/2 + O (e 2 ); in 
contrast, below a±, some of these will be attracted towards fixed point (i) where Asr — e/3. For these theories, the 
strong anisotropy exponent A will change dis continuously upon crossing the stability boundary between fixed points 
(i) and (ii). Note, however, that Asr — e/3 coincides with A^ = e/ (5 — 2a) at the line a = ax = e/ (3 + e) to 
all orders in e. Nevertheless, even though the exponents may undergo discontinuities, both the hybrid and the SR 
structure factors scale according to Eqn. (R). Thus, the scaling forms remain unchanged. 



IV. CONCLUSIONS 



Using field theoretic methods, we have analyzed the RG flow for a model of biased diffusion subject to a noise 
term, parameterized by its momentum dependence, q^ 1 a \ with spatially long-ranged correlations. One limit, a = 0, 
corresponds to a "short-range" model with purely local, conserved noise, and the other, a = 1, models biased diffusion 
with a non-conserved noise. The crossover from a < 1 to a = 1 presents no difficulties; however, the opposite limit, 
a — ► 0, is quite subtle. Here, the full crossover is observed only if we interpose a hybrid model with a = O (e), between 
the SR (a — 0) and the LR (a finite) theories. This hybrid contains the key elements of both, SR and LR, Langevin 
equations. It possesses a number of fixed points, including the equivalents of the SR and the LR models, so that the 
crossover can be understood in terms of a stability exchange between different fixed points. Clearly, considering a on 
the scale of e enables us to resolve this "fine structure" . Just below two dimensions, the scenario is following: for a 
below a lower stability limit ax, the theory is controlled by two stable SR fixed points, one FDT restoring and the 
other FDT violating, each with its own basin of attraction. Stated differently, we are in a "weak noise" regime, for 
d < d^R = 2 — 3a/ (1 — a), where the long-ranged noise does not significantly modify the universal behavior of the 
system. As a increases beyond ax but remains below an upper stability limit a 2 , the SR FDT restoring fixed point 
becomes unstable, leaving the SR FDT violating fixed point in control of the flow. This stability exchange is mediated 
by a fixed point line. Finally, above a 2 , the SR FDT violating fixed point also destabilizes and the nontrivial LR 
fixed point becomes globally stable. We note, in conclusion, that a similar mechanism, namely stability exchange 
through a fixed point line, has previously been observed in the Sine-Gordon model: there, two fixed points are stable 
below d = 2, corresponding to the high- and low-temperature phases, respectively. As d increases beyond 2, the 
high-temperature fixed point loses its stability, via a similar fixed point line pl| . 

Unfortunately, most of the nontrivial crossover phenomena discussed here are confined to dimensions 1 < d < 2. 
Above d — 2, the long-range noise dominates, either via its nontrivial or its Gaussian fixed point. In one dimension, 
on the other hand, there is no transverse subspace so that the SR FDT violating fixed point cannot be accessed. 
Consequently, the scenario described above must change significantly. It is conceivable that only the lower stability 
limit might survive here or that the two stability limits merge, leaving us with a stability exchange between the FDT 
restoring SR and the LR fixed point. Since, at a = ax, Asr = Alr — e/3 to all orders, this may be a reasonable 
conjecture. Moreover, in d = 1, we have a\ — 1/4 and A = 1/3. These values agree with the corresponding results 
p3| , p^[ for the one-dimensional KPZ equation with correlated noise where only one stability limit is observed. 

Nevertheless, our analysis plays the role of a pilot study for a number of other interesting problems. Clearly, one 
might consider an interacting theory subject to an external bias JlQj ] and a long-range noise term Jl6[ . Here, d c = 5, 
so that physical dimensions are more accessible, and comparisons with Monte Carlo simulations can be made. Other 
questions of interest concern noise terms with different spatial correlations, or nontrivial correlations in time. For 
equilibrium systems, the existence of an underlying Hamiltonian ensures that these correlations have no effect on 







static properties. For non-equilibrium steady states, however, the static behavior is gcncrically inseparable from the 
dynamics. Studies of anomalous noise correlations in non-equilibrium systems may help to unravel the nature of this 
coupling. 
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